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1 Introduction 



One of the main branches of theoretical and mathematical physics is the theory of exactly solvable 
models. The most successful approach to construct integrable two-dimensional lattice models of statistical 
mechanics and integrable (1 + l)-dimensional quantum field theory is through the solution of the Yang- 
Baxter equation Given a solution of this equation, depending on a spectral parameter A, one 
can define the local Boltzmann weights of a commuting family of transfer matrix r(A) |31 0] and the 
factorizable 5- matrix in a two-dimensional quantum field theory 0. 

The structure of the solutions of the Yang-Baxter equation based on simple Lie algebras is by now 
fairly well understood |6j. In particular, explicit expressions for the 7^-matrices related to non-exceptional 
afRne Lie algebras were exhibited in jT, and jSj. Since then, many other 7?.-matrices associated to higher 
dimensional representations of these algebras have also been determined j^]- 

A complete understanding of the vertex models living in a planar lattice include the exact diagonal- 
ization of the row-to-row transfer matrices which can provide informations about the on-shell physical 
properties such as free-energy thermodynamics and quasi-particle excitation behavior. The most efficient 
method for achieving this is the algebraic Bethe ansatz HIIJ, though its coordinate version is usually 
more efficient for finding the energy spectrum of concrete models 11 1|. A long-standing open problem is 
the diagonalization of transfer matrices of vertex models associated with solutions of the Yang-Baxter 
equation. 

One possible method of finding the eigenvalues of a given transfer matrix is the so-called analytical 
Bethe ansatz 112] . This technique relies on the unitarity, crossing and analyticity properties of the transfer 
matrix and, in some cases, an extra amount of phenomenological input is also required. This method has 
been applied to some of the models which we are going to consider in this paper, more precisely for the 
systems Iin \ Gn \ oi^^ ^J^J. Unfortunately, the explicit construction of eigenvectors of the transfer 
matrix is beyond the scope of the analytical Bethe ansatz. The construction of exact eigenvectors is 
certainly an important step in the program of solving integrable models |15| 

The importance of the algebraic Bethe ansatz does not rely on the calculation of the energy spectrum 
of a given model, but to also to supply information on the nature of the eigenfunctions. Thus is crucial 
in the investigation of off-shell properties such as correlators of physical operators ^ as well as for the 
calculation of form- factors ^Hl ■ 

A unified formulation of the quantum inverse scattering method for lattice vertex models associated 
to non-exceptional Lie algebras has been developed in the last years by Martins and collaborators [171 
1181 119j . In their works the Yang-Baxter algebra is recast in terms of novel commutation relations among 

(2) 

creation, annihilation and diagonal fields. In particular, the solution of the twisted D„_j_]^ vertex models 
is accommodated in their unification by the solution of a sixteen- vertex model 1201 . 

In this work, we will describe a detailed account of this method which complements results on the 
literature ^2], by the investigation of the trigonometric vertex models associated with the affine Lie 
algebras Bn \ Cn \ Dn \ ^'2n ^^i^ ^2ri-i- The D^^i vertex models are out of the scope of this paper. 

The outline of the paper is as follows: 

In section 2 we present the models through their TZ matrices. In section 3 the eigenvalue problem 
of the transfer matrix is formulated. In section 4 we perform a detailed construction of the intertwining 
relation in order to derive the fundamental commutation relations. In sections 5 and 6 the eigenvalue 
problem is executed in full detail for one and two-particle Bethe states, respectively, emphasizing the 
subtleties of each case and developing the language used in the text. In sections 7 the multi-particle cases 
are solved. The section 8 is reserved for our conclusions. In the appendix the A^^^ vertex models are 
considered for sake of completeness. 
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2 The vertex models 



The search for integrable models through solutions of the Yang-Baxter equation has been performed by 
the quantum group approach in |22| . where the problem is reduced to a linear one. Indeed, 7?.-matrices 
corresponding to vector representations of all non-exceptional affine Lie algebras were determined in this 
way in 

Quantum 7^-matrices for the vertex models associated to the Bn\Cn \ Dn \ A2^j and A^2n-i affine 
Lie algebras as presented by Jimbo have the form 
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(2.1) 



where we have introduced a label I, ^ = 0, 1, n — 1 in order to work with the nesting structure presents 
in the nested Bethe Ansatz construction. For a given value of n, the label I identifies a particular 
vertex model among those models with n — I < n. We can name the label Z = as the ground and 
the remained ones as the layers in the nest build. Here, Eij denotes the elementary Ni by Ni matrices 



{{E^0)ab = S^aSpb), where TV, = 2{n - I) for c|^'2„ d|^^2, and A'^^' and Ni = 2{n - I) + 1 for B^^Z 



(2) 



and A, 



(2) 

2{n-l)- 



(n-l)- 



The Boltzmann weights with functional dependence on the spectral parameter A are given by 



4'^ (A) = (c^-g')(c^-eo, 

4'^ (A) = -{n'-i){o^-ii), 



and 
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where q — e^^** denotes an arbitrary parameter and a' — Ni + 1 — a. The sign functions Eq = 1 (1 < 



a <n ~ I), Cq, = — 1 (n — Z + 1 < a < 2{n — I)) for c|j^2; and = 1 for the remaining cases. 



Here ^ and a are given respectively by 
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in the remaining cases. 

These 72.-matrices are regular satisfying PT-symmetry and unitarity: 

7^W(0) = x^\q)V('\ t41{\) = V^MlW^l ^S(A)7^^'l^(-A) = x^\\)x^\-\), (2.7) 
where V is the permutation matrix: P^'^ \a) <8) = <8) la). 

3 The eigenvalue problem 

In the context of two-dimensional classical statistical systems, each Yang-Baxter solution 7^(A) is inter- 
preted as local Boltzmann weights of an integrablc vertex model on a square lattice of size L x L. A 
physical state on this lattice is defined by the assignment of a state variable to each lattice edge. If one 
takes the horizontal direction as space and the vertical one as time, the transfer matrix tl{X) plays the 
role of a discrete evolution operator acting on the Hilbcrt space Ti^-'"^ spanned by the row states which 
are defined by the set of vertical link variables on the same row. Thus, the transfer matrix elements can 
be understood as the transition probability of the one row state to project on the consecutive one after 
a unit of time. 

The main problem is the diagonalization of the tl (A) matrix for these lattice systems. To do this we 
request the algebraic Bethe ansatz where the row-to-row transfer matrix can be constructed from a local 
vertex operator Cai{\)^ the Lax operator, which is viewed as a matrix on the iV-dimcnsional auxiliary 
space Va, corresponding in the vertex model to the space of states of the horizontal degrees of freedom. 
Its matrix elements are operators on the L-product Hilbert space Ti*^^-* = V^^ , where Vi corresponds to 
the space of vertical degrees of freedom and i denotes the sites of the one-dimensional lattice. An ordered 
product of Lax operators defines the N"^^ by N"^^ monodromy matrix 

r(A) = CaL{\)CaL-l{\) ■ ■ ' Cal{\). (3.1) 

which can be written as an A'' by A'' matrix with entries 

Tii (A) = ^I'k! (A) ® (A) ® • • • ® 41^, (A) (3.2) 

fci,...,fci_i=i 

where c\^\x) are N hy N matrices acting on the quantum space Vn- 

The transfer matrix of the vertex model with periodic boundary conditions can be written as a trace 
of the monodromy matrix on the auxiliary space Va 

N 

tl{X) = Tr,[T(A)] = ^Tii(A) (3.3) 

i=l 

and the eigenvalue problem is defined by 

TL (A)* A(A)* (3.4) 

where the eigenfunction 4" is obtained from the action of the non-diagonal matrix elements of T(A) on a 
reference state. 

Here one uses the fact that the Yang-Baxter equation 

7li2(A - M)^l3(A)7^23(M) = 7^23(M)7^13(A)7^12(A - m). (3.5) 
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can be recast in the form of commutation relations for the matrix elements of the monodromy matrix 
which play the role of creation and annihilation operators. The commutation relations are derived from 
the global intertwining relation 



^(A - ^l)T{X) (g) Tin) = T{p) ® T{\)S[\ - ii) 



(3.6) 



where we have used that the intertwining matrix S{X — fi) is defined on the tensor product Va (!S)Va and 
satisfy the relation S{X — ^) — VTZ{X — fj.) when the auxiliary space Va and the quantum space Vi are 
equivalent and the Lax operator identified with the TZ matrix, i.e., C{X) = TZ{X). The indices in the 
matrix TZ denote the spaces where its action is not trivial. 

In this paper a sufficiently general recipe is supplied to derive the fundamental commutation relations 
among the monodromy elements for the trigonometric vertex models associated with the Bn \ Gn \ \ 

(2) (2) 

and A2„_i afiine Lie algebras. 
For the vertex models H2.1(l the monodromy matrix (|3.3|l can be written as an Ni by Ni matrix 
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The usual reference state 



ioL)^'^=n®io), 



|0). 







V / 
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where Ni is the length of the vectors jO),-, is the highest vector of r(')(A) 

4"(A)|0l)'" = x['\x)\Ol)^'\, 

B^^HX)\QLf^ ^ {0, |0l)(')}, 

a ^ P=l,2,...,Ni 



dS(A)|Ol)(') = o, 

C«(A)|Ol)("=0, 



where 



x«(A) = [x^(A)]^ xr(A) 



(0/ 



[x«(A)]^ 



This triangular property suggests to write T'^'^(A) as a 3 by 3 matrix: 



4\X) Z?W(A) Bn,-iW 
C(')(A) P(')(A) B*(')(A) 



where one can identify four scalars 



(3.7) 



(3.8) 



(3.9) 
(3.10) 

(3.11) 



(3.12) 
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as well as, four {Ni — 2)-dinicnsional vectors 



B 



W(A) = (i?f>(A),sf (A),... , 13<'\X) 



I <!(A) \ 



(3.13) 



C(')(A) = 



(A) 



besides an {Ni — 2 ) by {Ni — 2) matrix denoted by 



P«(A) = 



.N,- 



D 



(0 



-2(A) / 



(3.14) 



(3.15) 



The problem of diagonalization of the transfer matrix becomes 



.(0 



Ni-2 



(0 - a(0 



Ai!'(A|{Aa)*«. 



(3.16) 



where ^I'm = (Ai, Am) is a scalar function named the m-particle state for the eigenvalue problem. 
In particular, the eigenvalue of the reference state (|3.8() is 



Af (A|0) = (A)]^ + (TV, - 2)[4'^(A)]^ + [y];;^.,(A)] 



(3.17) 



In order to construct other eigenvalues, one has to find the commutation relations among the elements 

of r(')(A). 
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4 Fundamental Commutation Relations 



In contrast to what happens to the six-vertex model and its multi-states generalizations it is a rather 
complicated task to find the commutation relations among the matrix elements of the monodromy matrix 
in a general case. However, as we are going to see, this construction can be simplified because all 
informations about the commutation relations are already encoded in the simplest cases. 

For Ni > 3 the corresponding Ni by Ni intertwining 5^'^ matrix can also be suitably written as a 9 
by 9 matrix in the form 
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VNiI 
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(4.1) 



where one can identify four iVj'^j^ -dimensional vectors 
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1=1 



(4.2) 
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N,-2 



(4.3) 



and a 7Vf_^^ by iVf.^^ matrix f(') = 'p{i+i)Y(.'\ with F^') obtained from the 7^(') matrix (EU by the 
reduction 
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(4.4) 



The remaining entries of [S"*-''] are five scalars {a;^'"' , yj']^ , yj;'^^ , yj^^i , J/^^^at, } and three A^j+i by Ni^i identity 

matrices , , 0:4 "* }. The i?i in the definition relations of the Y^''^ are column vectors of length Ni+i 
with only unitary element at ith position and t means transposition. 

Here we note that although these vectors and matrices act on the layer I + 1 (A^z+i = Ni — 2), their 
Boltzmann weights are written in term of the l-th layer Boltzmann weights. Therefore, our label I in a 
particular Boltzmann weight is indicating the model that it belongs. 

Now, using (|4.1|l and (|3.11|) we can write the fundamental relation in the form 



[^(')](A - fi)T^'\\) (g> t^'Hm) = t'^^Hp) ® r^''(A)[^''^](A - ^l) 



(4.5) 



in order to get 81 equations for the commutation relations among the entries of the monodromy matrix. 
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Before we look at these equations we would like to show what happens in the simplest cases. 
Ni — 3, all entries of [5*^'^] are reduced to the scalar status and we have 
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(4.6) 



This is the intertwining S matrix for the B^^^ and vertex models whose 7^-matrices are given by 
1)2. l|l with iV = 3. Moreover, the matrix elements of (|3.11|) are scalars in the auxiliary space 



(2) 



Ai Bi B2 
T = \ Ci B3 
C2 C3 A3 



(4.7) 



Substituting (|4.(j|l and H4.7|l into (|4.5|l we get the commutation rules for the entries of (|4.7|) proposed 
previously by Tarasov in the context of the Izergin-Korepin vertex model [23) . 

It is also worth note that the case Ni = 2 can be added to our discussion. In this case all entries 
with tensor status are removed from [S^''^. The result is 
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Xl 
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which is the intertwining S matrix for the C\ ,D 



and ^ 



vertex models whose 7^-matrices are also 



given by ()2.1)l . Consequently, their monodromy matrices preserve only the scalar entries ()3.12)l : 



T 



Ai Bi 
Ci ^3 



(4.9) 



Of course this is not a simple limit from the general case because their Bethe states are different by 
construction. However, we can derive the commutation relations for the entries of ()4.9)l by vanishing all 
entries which are vector or matrix in the general commutation relations which we are going to derive. 
Note also that C^^"* and A^^^ are six-vertex models while D^^'' is neither a six-vertex model nor a regular 
model as we can see from they TZ. matrices. 

In order to derive these so important commutation relations we shall proceed in the following way: 
we denote by E[i^j] = the (i,j) component of the matrix equation ()4.5)> and collecting them in 25 
blocks B[i,i], [i = 1, 5, j = i, 10 — i), defined by 



= E[10 - J, 10 - z] 



(4.10) 



For a given block B[i,j], the equation can be read from the equation Eij ( and the equation Eij from 
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the equation ) by the interchanging 

C(')(A) ^ C*^'H^i), Cn,-iW ^ CN,-iif^), V^'\\)^V^'\fi), 

^ ^3 ' yi,Af, ^ VNuI -"^21 ^ -^^12 ' ^ Ni,2 ^ " 2,Ni \^-^^) 

Here we note that the Boltzmann weights (|2.3|) satisfy the relation J/^'''^(A) — Up) ^'W- Taking into 
account these identifications, the computation to find commutation relations is considerably simplified. 
For instance, in the block B[l, 4] we can solve the equations En ~ and ei4 = in order to find 

A«(A)B(')(^) = z^'Ht^ - A)i?«(M)4'^(A) - 4|^^S(0(A)4')(^) (4.12) 

and 

4\x)B<^Hp) = z(')(A - m)B*«(a.)4'^(A) - "4|(A-m) ^.(0(a)a(')(^). (4.13) 

4^(A-/i) 

where we have introduced the function 

since it will appears many times in the text. 

Let us consider here two further blocks: the equations in the block B[2, 7] yield the following inter- 
twining relations 

vW 



4')(A)6*«(/i) - ^i^i^6*(0(^)40(;,)_^(0(A)®2?(0(^)J^^X>_^ 

^;jv,(m-a) 



B^,_^(^)C(0(A)-^lf>^i?^,_i(A)CW(M) (4.15) 
2/Ar,Ar, (/^ - A) ?;]v;^^(Ai-A) 



and 



4'^(A)i3(')(M) = f ^"-^\ g^')(^)4'^(A)- 5^^" ^^ g*(')(A)^p(')(M) 



+ 4^^^^Sa^,-i(a^)C*«(A) - 4?i^^S.,-r(A)C*(')(^) (4.16) 

where the tensor structure coming from (|4.5|) with scalars, vectors and matrices defined previously. It is 
also worth note the correspondence between H4.15|l and (|4.16f) via the exchange property H4.11|l and the 
matrices product order. 

A very important information is given by the commutation relations derived from the block B[2, 5] 
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4'^(A-/i) 



N, 



A) (A 



4'^ (A -A*) 



4?2(A-m) 



2/jv,w,(A"M) 

(4.17) 



and 



4 (a^-A) - A) y]Y,Ar,(A - m) 



(/i - A) ~ ' " ' x)^' {p. - A) TV, (a^ - A) 



-^§^^!-^B^\X)Af{^) (4.18) 
y]v w,(M- A) 



where we have defined a -/V;^]^ by f^l^^ matrix 



.W(A) = f(')(A)-^-^f«,(A)®4(^,(A) 



2^jv,jv, (A) 

which satisfies the permuted version of the Yang-Baxter equation 

42^(A)42(A + A^)42'(A) = 43'(A^)42^(A + A^)43^(/i) (4.20) 

In the definition H4.19|l the entries of s^'-' (A) were written in terms of the Boltzmann weight of the vertex 
model with label I. However, due to (|4.20() . its label must be A + 1. To write the matrix s(A) with its 
Boltzmann weights labeled correctly we can use the identities 

s(i)(A) 5('+i)(A) , , 

^ ' ^ ' a =1,2. (4.21) 



(A) 4'+^) (A)' 

where 5^'+^-' = p('+i)7?,('+^) and 7^''+^^ is given by H2.1|) replacing I hy I + 1. These relations give the 
emphasis of the meaning of the label I. Moreover, in (|4.17|l and (|4.18() we have used the following matrix 
identities 

^P^/HX-,) - -n^^(A-.) + #^y^?.(A-M) 
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and the scalar relation 

x<'>(A) 4"(-A) 

Many other commutation relations will be used in this paper. In the right time we will derive them 
recalling the block B[i,j] once more. 



5 The one-particle Bethe state 

In the quantum inverse scattering method the eigenstates of the transfer matrix are constructing by 
the action of the creators operators on the reference state. Such procedure results in excitations with 
multi-particle structure, characterized by a set of rapidities {A^} which are determined by solving the 
Bethe equations. In the general case the vector has Ni — 2 components and is used to define the 

one-particle state which is a scalar function obtained by the linear combination 

^«(Ai) = S(')(Ai).Ff^ Yl) (5.1) 

where is a vector with Ni^i components 

The action of the transfer matrix rf^^ (A) on this state is determined by (|3.9|) and the intertwining 
relations H4.5|) . The components of (|4.5|) needed for the construction of the nested Bethe ansatz of the 
one-particle state are the commutation relations H4.12|l . H4.1(j|l and (|4.17|l . In particular, the matrix 
relation H4.17() must be written in terms of its components in order to get the commutation relation for 
the scalar Daa{X): 

i?«(A)e«(^) = g(')(M)Tr.[ ^;'/^^^ ; y '''w]- "tn, ^y ^m^Ht^) 



a=l 



-i|%^[s*(0(A)«i(')]4')(^) 



where 1*^'^ is the A^/+i by Ni+i matrix identity and Tra is the trace in the auxiliary space. In 15.2|l . we 
also have used H4.21(l to write the matrix s^'^ as 5'*''+^-' and identified its permuted matrix 7^''+^^ with 
the Lax operator . 

In the nested Bethe ansatz procedure we always begin at the ground level / — and from now we 
shall omitte such label in the expressions of the first eigenvalue problem. 
The eigenvalue problem is accomplished by the action of Ti(A) on ^'i(Ai) 

tl(A)*i(Ai) - (^Ai{X) + i?„a(A) + A3(A)j S(Ai)^i |0l) 

= Al(A|{Ai})^i(Ai) (5.3) 

Taking into account the commutation relations H4.12|l . (|4.16f) , H5.2|l and H3.9|l we are able to turn the 
operators Ai{X), Daa{X) and ^3(A) over the creation operator B{Xi) and as result we have the following 



^'\X-^) yS(A-A*) 
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expression for the first eigenvalue problem: 



rL(A)vI/i(Ai) = Xi(A)z(Ai-A)v|/i(Ai)+X3(A)^^^^-^*i(Ai) 

yNN[>^ - Ai) 



-B(Ai)Tr„ 



4^(A-Ai) 



^i|0 



L 



-S(A)[X,(A,)^^4^ +X.(A,)^;4^m 10.) 
a;2(Ai-Aj a:;2(A-Ai) 

_ yAr2(A-Ai) ® l[X2(Ai) - Xi(Ai)]^i |0i) (5.4) 

yNN(A - Ai) 

Of course, the terms proportional to ^'i(Ai) are the wanted terms and contribute to the eigenvalue 
Al(A|{Ai}). The remaining ones are the so-called unwanted terms and they have to be eliminated by 
imposing restrictions on the rapidity Ai. 

In this case it is easy to see that we have two unwanted terms which are directly eliminated by the 
condition Xi{Xi) = X2(Ai). However, the trace term on the right hand side of (|5.4II does not give its 
wanted part directly. 

By simple inspection of the 7?.-matrices (12.1(1 for the ci^^D^^'' and Ai^rl-i vertex models, one can see 
that they trace in the auxiliary space is proportional to the N hy N identity matrix I 

Tr44V(A)] = (a:«(A) + (iVi-2)x«(A)+j;^^V^(A))l, iVi - 4, 6, 8, ... (5.5) 

It means that \I'i(Ai) is the eigenstate of tl(A) with eigenvalue 

Ai(A|{Ai}) = Xi(A)z(Ai -A) + X3(A)-'''^^"^'^ 




provided that 



-X2(A) (^-i^ + (7Vi-2) + -f^ ) (5.6) 
^i(Ai) 



= 1 (5-7) 



where we have used the notation 

^i'^(A)=2/^W(^-^i). XW(A)=.T«(A-Ai), a =1,2 (5.8) 

Here we make some remarks in respect to (|5.6f) . Although its computation is made on the ground (I — 0), 
the result also contains Boltzmann weights of the model in the first layer {I = 1). Of course, we can 
recall H4.21|l to write H5.6|) with all Boltzmann weights of the model on the ground. Nevertheless, as we 
will see later, this layer formation is very important -even when there is no nest to build. For instance, 
the eigenvalue (|5.6|) is not valid for the D2^-' vertex model because its first layer I — 1 (the d'^-* model) 
is not regular. It is also curious to note that no further constraint is necessary for the vector J-i. 

However, when Ni is an odd number, the situation is not so simple because the Tra[-C^''] is not more 
proportional to the identity matrix due to the weight 2/^'^(A) with a ~ P and a — a' (|2.3|l . the element 

common to the diagonals of 7?,*^'^. 

The trace in ((5.4|l represents the transfer matrix of i + 1 site chain with one inhomogeneous site 
coming from the Lax operator (A — Ai). To solve the eigenvalue problem of r|^j(A, {Xi})Ti \0l), we 
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first note that t^I^ G H^^^ ^H^'^'^ and the part of t|^|^ hivolving 2?(A) G H^^^ commutes with Ti G H^'^^ 
and can hit directly the reference state |0l) 



,(A, {Aiim \0l) = X,iX) (rf )(A, {Ai})^i) |0i) . (5.9) 



It means that the eigenvalue condition for (|5.4() leads to the requirement that J-i ought to be an 
eigenvavector of t'^^\x, {Ai}). Therefore, if we suppose that 

)(A, {Aiim = a['\x, {Ai}| • • • (5.10) 
the last wanted term in (|5.4(l has the form 

X,{x /''^^J^^;^l--^ B{X,)T, 10.) (5.11) 

and we conclude for Bn'' and A'^^ vertex models that ^'i(Ai) is an eigenfunction of tl{X) with eigenvalue 

A.(A|{Ai}) = Xi(A)z(Ai - A) + X2(A)^^ii^iMAl^ + Xs{X)^^^^-^ (5.12) 

X^^'{X) yA'Jv(A-Ai) 

provided that 

MxT) = ' ^'-^'^ 

The eigenvalue H5.12|l is partial because we still have to solve the eigenvalue problem H5.10(l in order to 
know the value of A^^-'(A, {Ai}| • • • ). Here we have reached a point which is typical of nested Bethe ansatz 
problems. It means that we have to solve an another eigenvalue problem for the transfer matrix t[^^ 
with its 7?,-matrix given by H2.1|) but with ^ = 1, i.e., the first layer of the nest. 
The row-to-row 1-site inhomogeneous transfer matrix Ti^\x, {Ai}) is given by 

Ni-2 

4'\X) = a['\x) + i?«(A) +^«(A) (5.14) 

The notation used in (|5.14|l is to be considered as a shorthand as these terms depend also on the inho- 
mogeneous parameter Ai. The reference state jOi)^^-* is given by (|5.9|l and it is a highest vector 

A['\x)\0,f^ = xW(A)||Oi)W, i?W(A)|Oi)(^)=xW(A)|Oi)«, 
4^)(A)|0i)« = X«(A)|Oi)«, Z?S(A)|Oi)«^0, 

c!^H\)\Oif^ = 0, i?«(A)|0i)(^V{0,|0i)(^)}, a^P = l,2,...,N,-2. (5.15) 
with eigenvalue 

A«(A, {Ai}|0) = X['\X) + im - 2)X^'\X) + X^'\X) (5.16) 

where Xl^^ (A), a = 1,2,3. are given by 1)5. 8|l . 

The condition that Ti ought to be an eigenvector of t[^\x) requires the diagonalization of t[^\X), 
which can be carried out by a second Bethe ansatz. There are two candidates for Ti G 7Y^^^: the own 
reference state |0i)''^^ and the one-particle excitation S^^^(A^^^) |0i)'"^'. 

The choice !Fi = |0i)*-^'' implies in a particular linear combination for the one-particle state (|5.1ll and 
in this case \l/i(Ai) is an eigenfunction of tl(X) for the B^^^ and A2^'_^ with the eigenvalue expression 
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equal to H5.6|l and the Bethe equation equal to H5.7|l . However, the second choice seems to be the most 
general but we will not consider it now. This case will be presented later in a more general context where 
the multi-particle state is treated. 

Notice that we have not yet given an explicit rule to eliminate the unwanted terms, until now they 
were merely canceled out. In order to find such rule we will consider the two-particle state with more 
details. 



6 The two-particle Bethe state 

In analogy with the scalar case |23[I24| we have two types of linearly independent contributions yB*^'' (X'/B^'-' 
and -Bat,-! for the vector $2 ''(Ai, A2). We seek vectors in the form 

$^'^(Ai, A2) = B^'\Xi) e(''(A2) + Bw,„i(Ai)r(Ai, A2) (6.1) 

where r(Ai, A2) is an operator-valued vector which has to be fixed such that <I'2^(Ai, A2) is unique. 

It was demonstrated in |23j that ^2\^ii ^^2) is unique provided it is ordered in a normal way: in gen- 
eral, the operator- valued vector $m (Ai, . . . , Am) is composite of normal-ordered monomials. A monomial 
is normally ordered if in it all elements K*^'-' , and Bn,-i are on the left, and all elements C''-', C*^'^ 
and Ctv,-! are on the right of all elements A''^\ V^^^ and ^3''. Moreover, all elements of one given type 

have standard ordering, i.e., T^?^, (A^T^'^^^JAa) • • • T^'l/3„ (^'»)- 

Now we recall the intertwining relation i|4.5() to solve the equations of the block i3[l, 5] in order to get 
the following commutation relation 

+5A.,_,(A)A(')(^)^Mi^, (6.2) 

from which we can see that (jB.ljl will be normally ordered if only if it satisfies the property 

$«(A,,A2)=<i>('^(A2,A,)^Ji^5J^^i-^ (6.3) 

2;r''(Ai - A2) 

This condition fixes r(Ai, A2) and our vector for the two-particle state has the form 

.(0 



1^ (\ — A ) 

(Ai, A2) = 6(')(Ai) ® SW(A2) + Bn,^^{X,)4\X2)^^ (6.4) 



Here we notice that the condition (|6.3|) must be generalized to include multi-particle state and it will 
play the main role in the elimination rules of the unwanted terms. 

The action of the transfer matrix on this vector is more laborious. In addition to (|4.12() . (|4.15l) and 
H5.2|l we appeal to (|4.5() to derive nine further commutation relations. 

Due to the presence of the scalar Bni-i in H6.4|l we have solve the block B[\, 7] of 14.5|l in order to 
get its commutation relations with Ai and A3 in a given layer I 

4\X)B,,-1{,) = ^%^BA.,-i(/.)A«(A)-4f^^BA.,_,(A)4'\A*) 

VNMii^-x) 4;j^^(/i-A) 

-gW(A)0gW(A.) ^f'^^~^\ (6.5) 
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and 



4"(A)i?A.,_i(/i) 



s*("(a)®s*(')(m), 



(6.6) 



The block B[2, 6], i3[4, 6] and i3[4, 8] can be used to derive the relation 



(6.7) 



Again, we must work out H6.7I in order to get the commutation relations of i?Ar,-i with the scalar 
Tr,[pW(A)] 



Tia iX)]BN,-i{fi) = (A - ^i)z^''^ {^i - A)S^-i (/^)Tr, [p(') (A)] + l^BN-i{X)TTa [V^'^ 



4^(A-/i)2 



xf{X-ti) 



x'i>{x-^i) 



(^(0(A)i3*(0(^) _ Tr,[6*(')(A)B(')(Ai)]) 



where we have used the identity 



1 - 4^i>^-^^^'(>^-^^ ^ ,(0(A - ^),(0(^ _ A) 



x^\x-^Ji) xf{x-i,) 



(6.8) 



(6.9) 



Note also that these scalar relations will survive the reduction mechanism for the six-vertex models 
previously presented. 

Since we need a second commutation step in order to that t|^'' hits the reference state, we compute 
the commutation relations among creation and annihilation operators 



C«(A) ® S«(/z) = B«(/z) ® C(')(A) - ""f,}^ 

x('^(A-m) 



A^l^ (A)-D(" {h) - (/z)2)(') (A))] (6.10) 



C«(A)i3,,.,(;.) = 4^i^S.,„,(^)C(0(A) + 4^i^^B*(')(A)4')(^) 

y['i{fi~X) g*(0(^,)40(;,)_p(0^g(0(A)C^,_,(^) ^(6.11) 



both relations were obtained from the blocks B[2, 2] and B[4, 7], respectively. We also need the commu- 
tation of i?Af,_i(A) with the other creation operators 

i3^,_i(A)BW(/.) = z(')(m - A)B(')(Ai)B^,_i(A) - " ^) g(0(A)i?^^„,(^), (6.12) 

4'(m-a) 
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(/) 

B*W(A)B;v,-i(m) = z^^\X - M)Siv,-i(M)B*«(A) - ^±^^^Br,,-i{X)B<^\i,). (6.13) 

which are obtained from the blocks B[l, 6], B[l, 8] and B[l, 10]. 

Here we observe that the final action of t^^ (A) on normally ordered vectors must be normal ordered. 
This implies an increasing use of commutation relations needed for the eigenvalue problem. For instance, 
the action of the scalar Af\X) on the vector $2''('^i) ^^2) has its normal ordered form given by 



4'^(A)$^'^(Ai,A2) = n^W(Afe-A)$^')(Ai,A2)4'^(A) 



fe=i 



4^(Ai-A) 

4^(A2-A) x\ '{X1-X2) 

+Sjv,-i {X)G^^I (A, Ai , A2)4'^ (Ai)4'^ (A2) + • • • (6.14) 
where we have defined a row vector with N^^-^^ entries 

G«(A Ai A2) = - X) Yj:! ,{X - X,) s(^+^HXi ^ X) ^ y[%Mi-X) y}:^ .{X, - X,) 

' ' 4'^(A2-A)y(J)^^(A-Ai)4'+^)(Ai-A) y«^^(Ai - A) 2/«^,(Ai - A2) 
which satisfies the cyclic permutation property 

gS(A, Ai, A2) = g(')(A, A2, Ai)^fJJJJi^l^ (6.16) 

a;r'^(Ai-A2) 

Observe the mix of layers with respect the Boltzmann weights of G21 which really acts on the I + 1 layer. 
The action of the scalar A^\x) on the vector $2^(Ai, A2) has a very similar normal ordered form 

aW(A)$(')(A„A2) = n f ^)~':\ ^^'^(Ar,A2)4'HA) 

k=i yjv,Jv,(A ~ Xk) 

-^(')(A2i) ^^'^^,^~y, B*«(A) ® BW(Ai) ® 1?W(A2) 

VNiNliX — X2) 

_,(0(A,2)^|^i^^ s*«(a) ® eW(A2) ^ p«(Ao4£Jr^^^ 

y^U(A-Ai) 4'+^^(Ai-A2) 
+B;v,-i(A)h('i)(A, Ai, A2)©W(Ai) ® 2?«(A2) + ■ • • (6.17) 

where the vector H21 is given by 



hWqa A^ yNgi^-^^) yff2(Ai-A2) x^')(A-AO I^^Ua - A2) 

H2i(A,Ai,A2j - — — — - — — - (,) — rr (0 — TT i^-i^) 

yAr,Ar,(A - Ai) y^, jv,(Ai - A2) VninM ~ ^1) 2^JV,JV,(A - A2) 
and satisfies the cyclic permutation property 



H<\^(A,Ai,A2) = hW(A,A2,A0 '^!,'|]!|^^ (6.19) 

xr'^(Ai - A2) 
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Finally, the action of the scalar TTa['D^^\X)] on the vector <i>2''('^ij ^2) is a little bit different 



Tr. [P(') ( A)] 0,(0 ( A, , A. ) = ( A„ A. )Tr„ { ^[^1^' 4^77^^'^ ^^^^ 

X2 (A — A2) X2 (A — Ai) 

_,(0(A, - A.)4J^^^bW(a) 8 S«(A2)[P»(A,) « i(0] 7^<'+^HAi-A.) 

4^(A-Ai) 4'+^^(Ai-A2) 

_,(0(A, - Ai)4r^^i3(0(A) ® S«(A0[p(0(A2) ® iW] 



+.("(A: 



4'^(A-A2) 
.(0 



+,(0(A, _ A2)^f^^^^[i?*W(A) ® S(0(A,) « i(0]^g^J(^1^40(A,) 
yj;U(A-A2) xf+^)(Ai-A2) 

+i?Ar,_i(A)Y^'i)(A, Ai, A2)[I?(0(A2) 1«]a(0(Ai) 

+i?^.,_l(A)Y(0(A, A2, Ai)[P('HAi) l«]4)(A,)^^^Ji^l^ + . . . (6.20) 

a;r'^(Ai - A2) 



where the vector Y2';^^ is given by 



Y(')(A, A„ A2) ^ [.(0(A - AO^fei^ - 

4''(A-A2) 4^(A-Ai)x«(Ai-A2) y!^^v,(A-Ai) 



(6.21) 



and also satisfies the cyclic permutation property 

Y^'>(A, Ai, A2) ^ Y<',\a, A2, Ai) ^|'|'|^^^ ~ (6.22) 

4+'^(Ai - A2) 



The ellipses in these expressions denote normally ordered terms containing factors of the type C'-O, C*(0 
and Cni-1- 

It should be worth note that we have used some matrix identities to derive (16 .141) - (16 .2011 



5('+i)(Aafc)5('+i)(Afc,) 



4^^^ (A.,) 4'^'Hx,a) 

4')(A,,)^('+l)(A.fc) y2M(Aab) Yj^liXac) ^ xi;\Xa,) xi;\X,a) ^ ^ xi;\X,,) S^' + ^H^c) 

x'^^H^b) X^2^^\Xab) ^^^^(Aab) y^Af(AQc) X^'^ (A^b) S^'' (Acq) 2;^'^ (Acb) X^'^^^ (Aac) 

3^2^ (Agfa) yif,2(Aac) , X^'^(Acfc) yif;2(Aafc) _ m,. ^ yjr^UAac) / J, / ^\ (-R OQA 

where / is a Ni_^_^ by Ni_^_^ matrix identity and Xab — — A^. 

Let us begin with the eigenvalue problem for the two-particle state in a homogeneous L site lattice 
which is defined by the linear combination 

*2(Ai,A2) = $2(Ai,A2)^2|0l) (6.24) 

where $2 is given by H6.4|l and the vector T2 has components f"^ G C, (a, /3 = 1, iV;+i). 



16 



The action of tl{X) on ^'2(^1, A2) can now be computed using Ht).14|l . (|6.17|) and H6.2U|I : 



2 



tl(A)*2(Ai,A2) = Xi(A)[]z(Afe-A)*2(Ai,A2)+X3(A) J]^^^:^-^ *2(Ai,A2) 

k=l k=l y^'^^ ^' 

+X2(A)$2(Ai,A2)^3^ [t!^\\)T2} 

~^4t4^W ® 'i'i(A2)[Xi(Ai)z(A2i) - X2(Ai)z(Ai2)^jSJ^^]^2 |0l) 
a;2(Aio) a;i^^(Ai2) 

_X3(A20)^ ^ $^(Ai)[Xi(A2)z(Ai2)4fr^ - X2(A2)^(A2l)]^2 |0l) 

a^2(A2o) a;i^^(Ai2) 

_ iA^Ao2) g,^^^ ^ ^^^^^^ ^ l[X2(A2)z(A2l) - Xi(A2)z(Ai2)^i^^]^2 |0l) 

yNN{Xo2) a;i^^(Ai2) 

_ iWAol) ^ ^ ^ i[X2(Ai)z(Ai2)4fr^ - ^i(Ai)z(A2i)]^2 |0l) 

yjviv(Aoi) a;i^^(Ai2) 

+BAr_i(A){G2i(A, Ai, A2)Xi(Ai)Xi(A2) + H2i(A, Ai, A2)X2(Ai)X2(A2) 

+Y21 (A, Ai , A2)X2 (A2)Xi (Ai ) + Y21 (A, A2 , Ai)X2 (Ai)Xi (A2) \f }^2 |0l) 

xl (A12) 

(6.25) 

where we have used the rapidity difference notation with Aq = A and the definitions 

2 

$i(A,) -6(A,), X^^\X)^l[x['\X-Xk) (i = l,2) (6.26) 

k=l 

Moreover, the unwanted terms were combined and we are presenting the two-site inhomogeneous transfer 
matrix for the first layer 

4'\X) = Tr„ [/:«(A - A2)/:iV(A - X,)] . (6.27) 
Before we see how the unwanted terms can be canceled, let us first to consider the eigenvalue problem 

TW(A)^2=AW(A,{AJ|{Af)})^2 (^-1,2), (6.28) 

where our choice for the vector J-2 is implicit 

^2 = (f>«(AW,AW)|02)('^ (6.29) 

and above $2^'' is given by H6.4|l and |02)'"'^' by (|3.8|) . From these results we have an incomplete expression 
for the eigenvalue 



2 



X2(A-Afc) , ^ ,,,A«(A,{AJ|{A«}) 



A,(A|{A„ A2}) = x,ix) n ziXk - A) + MX) n "7, + x^ix r^ (^-so) 

fe=i ^JiVNNiX-Xk) xi'\x) 

because A2'^-'(A, {Ai}|{A^^''}) is still unknown. 

At this point it is worth giving further information about the shorthand notation used in this text: 
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• tl(A) is a row-to-row homogeneous transfer matrix in a L site lattice and Am(A|{Ai}) is the eigen- 
value of tl{X) for the m-particle state with rapidities A,;, i = 1,2, ...m. Acting with Ti(A) on its 
reference state |0l) we have the factors 



Tm (X) — ''n'l'' (A, {A^' '^'}) is a row-to-row inhomogeneous transfer matrix in a m site lattice with 
inhomogeneous parameters A^-' — 1,2, ...,m. Its eigenvalue for the m-particle state is denoted 
by Am (A, {Af ~^^}|{Af •*}) where A-'^ are rapidities of the particles and xf~^^ are inhomogeneous 
parameters of the corresponding lattice. Acting with Tm (X) with / > 1 on its reference state \0m) 
we have the factors 

m 

x«(A) = xW(A,{Ar)})=n-i'(^-^r^)' «=i'2 

rn 



X,(A) = [xa(A)]^, a =1,2 



and 



^3 (A) = [vnn{\)Y 



(6.31) 



.32) 



k=l 



where m is the particle (site) number in the layer I and A^'^'' — A 



In respect to the unwanted terms of (|6.25() . the main motivation to write this section, we first recall 
the relation (|6.3|) in its generalized form by using the cyclic permutations of the factors in the normal 
ordered vector for the m-particle state in the l-th layer 



l-^l: ^2, A„J - ^)rl\^2, Am, -^ij ^ ^ — (1+1),^ TT (i + l) , ~~ (b.66) 

Xi (Ai — Am) 



(A1-A2) xi'+'^(Ai-A3) 



Now we define an operator M by 
M([+i)(A,{AJ) = Tr, 



S^j^ \Xl — A2) 5*^2^ ''(•^1 ^ -^3) Sijn — Am) 



-r^' (Ai - A2) (Ai - A3) (Ai - Am) J 



0+1) 



(6.34) 



and we remark that M is the normalized permutation of a m site inhomogeneous transfer matrix 

r^'^A, {A,}) = Tr, [c^^H^ Am)4™-\(A - Am-i) • • • cT\>^ - Ai)] (6.35) 
Thus, we can write the generalization of H6.3II in a more convenient form 

<i>«(Ai, A2, Am) = $^'(A2, Am, Ai)Af(^+i)(Ai, {A,}) (6.36) 

where j = 1, m. 

Using these expressions with / = 1 and m — 2, one can see that the unwanted term B{X) $i(A2) in 
can be eliminated by the condition 



.(1) 



z(A2 - Ai)Xi(Ai) - z(Ai - A2)X2(Ai) — 



(Ai) 



.^2=0 



2^1 (A1-A2), 



(6.37) 
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and from the fact that B{X) (g) $i(Ai) is a permutation of B{X) <l>i(A2) by the interchange Ai <^ A2, 
one can use (|6.36|) with I — 1 and m = 2, in order to get a second ehmination condition 



^2''(Ai, {A.}) I z(Ai - A2)Xi(A2) - z(A2 - Ai)X2(A2) 1 ^2 = (6.38) 



4'^(A2-Ai), 



Consequently, we have the foUowing eigenvalue problems 



r(^)(A.)^2 = x«(A. - A,)4^^^J-2 (a ^ 6 = 1,2) (6.39) 

which should be seen as a generalization of the Bethe equations since this restriction alone eliminate all 
unwanted terms. To see this, one can use the same statements in B*{X)(E)^i{Xi)<^l and S*(A)(g)<i>i(A2)® 1 
terms of (|6.25ll in order to get the same pair of equations (|6.39|) . Finally, substituting (|6.39|) into the 
-BAr-i(A) unwanted term of H6.25|l one can see that it is also canceled out. 

It is now clear that the eigenvalue 16.3()|l as well as the Bethe equations are written in terms of the 
eigenvalue of T2^^ (A) : 



XijXg) _ z(Aa-Afc) A^^HA„{A,}|{Af^}) 

^2(Aa) z{Xb^Xa) x['\Xa~Xb) 



(07^6=1,2) (6.40) 



where A^^^(Aa, {A.}|{Af ^}) is the residue of A^^^ (A, {A.}|{Af ^}) at A = A^. Consequently, we have to 
consider this second eigenvalue problem (|6.28|) in order to fix the results of the first eigenvalue problem. 
Indeed it is a simple task because all computation was already made and the result follows directly read 
from H6.25() with the following trivial modifications: 



Changing the vertex models: (B^^\cL'\D5.^\A^t\A^t^_ J^(bW i,CW i,dW i,A^^i_,),A^^i_,)_ J, 
putting the label / = 1 in the new Boltzmann weights. 

• Changing the lattice: L site homogeneous lattice {tl(X)) 2 -site inhomogeneous lattice (t^^^(A)). 

• Changing the two-particle state: $2(Ai,A2)JS \0l) ^ <S?i^\x'-^\ xi^^)T^'^^ |02)^^^ 

After these modifications we still have to consider a third eigenvalue problem by repeating this chang- 
ing procedure by increasing by one unity the value of the label I but keep the size of the last lattice. 
Next, we repeat the last procedure and so on. Of course, such a procedure must have an end. 

In the nested Bethe ansatz procedure the end is obtained in the layer for which the eigenvalue and 
the Bethe equations are fixed. It means that we are work with the models in the last layer. 

The value I = n — 1 defines the last layer as being the B^^^'' and A2^' vertex models {Ni odd) and as 
the Cj^'' jDj^"* .an Ag^"* vertex models (TV; even). Thus they must be considered separately: 

6.1 Bn^ and Ag^^ two-particle state 

We already mentioned that these models possess a limit via the reduction to scalar status. It means that 
in the last layer we are working with (|4.6|l and (|4.7|l with the Boltzmann weights of the B^^'' and A'^^ 
vertex models. The explicit expression for the eigenvalue problem of the L site homogeneous transfer 
matrix for these nineteen- vetex models was already presented in the reference |24| . However, here this 
expression will be presented by a reduction procedure from our general formulation. 
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Consequently the eigenvalue expression Ht).25|l is reduced to 



k=l k=l k=l y^^^ 

a;2(Aio) V t^(Ai2)/ 
X3(A2o) f ^^(^^)f(Ai2) _ x,iX,)z{X,^)] Bi(A)Bi(A2) IO2) 



a;2(A2o) V ^(Ai2) 

-^^4^ f^2(A2)z(A2i) - Xi(A2)4^ ) S3(A)Bi(Ai) IO2) 

y33(Ao2) V t^(Al2) 

y32(Aoi) ^^^(^^)£iAi2l _ x,iX,)z{X,,) ] B3(A)Bi(A2) IO2) 



y33(Aoi) V t^(Al2) 

+B2{X) {G2i(A, Ai, A2)Xi(Ai)Xi(A2) + i?2i(A, Ai, A2)X2(Ai)X2(A2) 

+ r2i(A,Ai,A2)X2(A2)Xi(Ai) + r2i(A,A2,Ai)X2(Ai)Xi(A2)— ^1 IO2) (6.41) 

t^(Al2) J 

where we have omitted the label I — n — 1. 

In H6.41|l we have used the function a;(A) which is the reduction of the objects S, TZ and Tra[£] 
present in (|6.25|) 

5(A) _ 7^(A) _ Tt,[C{X)] _ 1 f^^^^^^ _ 2;23(A)y32(A) 



a;(A) a;i(A) a;i(A) xi{X) xi{X) \ yssW 

w(A)w(-A) = 1. (6.42) 
Therefore, for the last layer the eigenvalue is 

A2(A, {Af- ^)}|{A,}) ^ X,iX) n ^(A. - A) + X,{X) f[ ^^^-^ + X2(A) fl (6.43) 

k=i Li 2^33(A - Afc) ^J-^ a;(A - A^) 

provided that 

Xl{Xa) z{Xa-Xb) /\ \ N /, 1 o AA\ 

Y t\ \ ^~F\ r-:^i^b-Xa), a^6=l,2 (6.44) 

The factors Xa are given by 

2 2 
^3(A) = n y33(A - Ar') ,Xa{X) = l[xa{X-Xr'), (a = l,2) (6.45) 

fc=l k=l 

where the inhomogeneity parameter {A"^^} make the link with the next to the last layer. 

Let us go back to the ground in order to write the full nest through of a sequence of terms where the 
models are explicitly identified: 



A,(A|{Ai,A2}) = Xi(A)nz(A.-A)+X3(A)n + X2(A)-^^'^^^' 

k=i 



(0) e (B«,A(^)) and (l)e(Bi^2i,A(^t^.) 



^J-^yNN{X-Xk) ^2'^ (A) 

(2) N 
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Ar^(A,{A«}|{Ar)}) 



''2 

(0 e (B«,Agi_,)) and (/ + 1) £ (B«,_i, A^^) 

I = l,2,...,77-2 



A^"-^^(A,{a|"-^^}|{a|"'^^}) ^ A<"-^)(A) 



2 



X("-^)(A) -jg. x^"'^^(A-A<^"'^^) 

^ z("-i)(A-Ar^^) 
"-y, .(A-Ai"-)) 
e (BW,Af) (6.46) 

Remember that the inhomogeneity parameters {Xf are imphcit in (|6.46|) through of the definition of 
Xi^\x), a = 1,2, 3. The function w(A - A^""^^) is given by 
The corresponding Bethe equations are 



XiiXa) _ z(A,-Ab)AW(Aa,{AJ|{Af^}) 



X2{Xa) z{Xt-Xa) x^^\Xa-X,) 

Af)(Ai')) ^ z(')(Ai')-Ai'))A^^)(Ai'\{Af)}|{Af+^)}) 
X(')(A«) z(0(Ai')-Ai')) xr^)(A«-Al')) 
/ = l,2,...,n-2 



a 7^ 6 1,2 



where 



^1 , '^1' ,> (A"-^)-Ai"-^)), a 7^ 6 =1,2 (6.47) 

x("-^^(Ai"-^)) 2("-i)(A["-^)-Ai"-')) 



Ar)(A«,{Af'}|{Ar)}) = res Ar^(A,{Af'}|{Ar^)}) 



.r)(o)xr^(AW - A«) n -^'+^)(Ar^ - Ai')) 



fc=i 

I = 0,l,...,n-2, 07^6=1,2. (6.48) 

It is curious that the solution of the eigenvalue problem of the L site homogeneous transfer matrix for 
the two-particle state is given in terms of the eigenvalue problems of the two site inhomogeneous transfer 
matrices for two-particle state. We remark the participation of n different models in the construction of 
the nested Bethe ansatz for the B^^-* and A'^J^ vertex models. 
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6.2 C„ , Dn and A''2n^i two-particle state 

For Ni = 4, 6, ... the last layer involves the C2^'',D2^'' and Ag^' vertex models for which their nests are not 
complete. Indeed, just more one layer is necessary in order to complete them. To do this we recall H4.8(l 
and (|4.9I) in order to include the C^^^D^"'^^ and A^^-* vertex models in our discussion. 

For the Ci and A\ ' vertex models, the L site homogeneous transfer matrix is the trace of 



and the reference state 



Ti(A) = Ai(A) + A3(A) (6.49) 



|0l)= ( I) (E>---(^(1 ] (6.50) 



is a highest vector of 14.9(1 



Ai{X)\Ol) = ^i(A)|Ol), A3(A)|0i) -X3(A)|0l) 
Ci(A)|Ol) = 0, Si(A)|0i)^{0,|0i)} 

Xi(A) = [xl(A)]^ X3(A) - [y22(A)]^ (6.51) 

The two-particle state is a new state without any relation with ((6.4(1 which can be defined by 

*2(Ai, A2) = 4'2(A2, Ai) = Bi(Ai)Bi(A2) (6.52) 

The action of tl(A) on this state can be computed using the commutation relations ((6. 5(1 - ((6. 7(1 in their 
reduced form 

r,(A)*2(Ar,A2) ^ f^i W fl ^^^^ + ^3(A) fl ^^^^ 
\^ ^Jl V22 ( Afc - A) 2/22 ( A - Afe ) ^ 

y22(Ai-A)V 2;22(A2-Ai) y22(Al-A2)/ 

yi2(A2 - A) ( ^^(,^) ^i(Ai-A2) ^ X3(A2)^^^^') Si(A)i?(AO 10.) 




y22(A2-A)V 2;22(Al-A2) 2/22(A2-Ai) 

where we have used the property of the definition 1(6.52(1 and the identity 

yi2(A) y2i(-A) 



(6.53) 



?/22(A) y22(-A) 



0. (6.54) 



The eigenvalue is 



provided that 



A,(A|Ar, A2) = X,{X) n ^^^^^ + ^3(A) fl ^^^^ (6-55) 

2/22 (Afc - A) y22(A - Afc) 

XljXq) Xi{\a ~ Ab) j/22(Afc - Aa) / i o c;fi> 

Y f\ \ ^ —r\ r\ TT' a/ 0=1,2. (6.56) 

A;3(Aa) 2:1 (Ab - Aa) 2/22 (Aa - Afc) 

We now turn to the diagonalization problem of the d'^^ vertex model. It turns out, however, that the 
Lax operator of this model can decomposed in terms of Lax operator for the six- vertex model associated 
with the A^^^ Lie algebra. It means that from the isomorphism D2^-' =A^^-'®A^^-', we can write 

C°'iX)=C^'{X)(g>C^'{X) (6.57) 
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Here, a more careful analysis is required. First we identify the Lax operator C^^ with the corresponding 
7^- matrix H2.1|l and then we make a sign transformation in the Boltzmann weights Uafii^) that preserves 
the spectrum of the transfer matrix associated i.e., yafii^) ^ ^2/a/3(A) for a ^ [3 and a ^ (3 . Now is 
not difficult to verify the isomorphism (|6.57|) where C^^ (A) is identified with the 7^-matrix of the A^^^ 
vertex model listed in the our appendix. 

Consequently, the eigenvalues of the model 02^-* are given in terms of the product of the eigenvalues 
of two A^^-* six-vertex models: 

AfMA|A±,A2±) = A+(A|A+,A+)A-(A|A^,A2-) (6.58) 

where 

2 2 

A±(A|Ai, A2) = Xi(A) n ^(A^ - A) + A2(A) [] ^(A - A±) (6.59) 

fc=l k=l 

provided that 

^iIA|2 ^ z^A|-A|2 ^^^,^12. (6.60) 

X2{Xt) z(A±-A±)' ^ ' ^ ' 

Here [X^ , A^} and {Aj~, Aj"} are rapidities of the two-particle state related to each one of the two six- 
vertex models. The algebraic Bethe ansatz for aI^'' vertex models is discussed in the section 8 of this 
paper. 

From these results we can see that the two particle state ^'2(Ai, A2) is an eigenstatc of the homogeneous 
transfer matrix Ai(A) with eigenvalue 

2 2 
Ai(A|{AJ) = Xi(A)nz(Afe-A) + X3(A)n^^^^-^ 

+X,{X) (£G(')(A,{Af-^)}|{Af^}) + r ) (6.61) 



where 



^2 i-^J fc=l ^2 W k=iyNiNii^ - \ ) 

and the last terms T depend on the model: 

T-G. (A.iA. }|(A. » + :^T^ li - A) + .<»>(A - A«) 

for the Cn^ and A2^^_j^ vertex models and 



T = [Aj"-^)(A)n-^"-^HA+-A)+Ar''(A)n-^"-'HA-A+)] 

k=l k=l 

X[X[-'\X) fl z("-i)(A,- - A) + 4"^'^(A) n ^'""'HA - A, )] 
fc=i fe=i 

(n-1) G A^^^ (6.64) 
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for Dn^ vertex model. The last row in (|6.64() is to remember that the Boltzmann weights are those 
presented in the appendix. 

The sequence of the Bcthe equations has the form 





XliXa) Z{\a~ 


Xb) 




X2iXa) ziXb - 


Xa) 




) z«(A«- 






) -('HAi'^- 


A«: 




I = l,2,...,n- 


- 2. 



i^)(0)n^''^(AL''-A„) a^6=l,2 
fc=i 

(0)n.('+i)(Ar)-A«), a^6=l,2 



(i+i 

fc=l 



(6.65) 



which will end in two different ways: 



A^tM^^t!!) - ^^-^HA.^-A,^) 



e A^^^ (6.66) 



for d1^^ and more two equations 



x("-i)(Ai"-'^) z(»-i)(Ai"-'' -Ai"-')) ^^11 



(n-f) e (C«,A(^) 



fe=i 

2 ' ^3 



( Ai") ) _ ( Ai") - a[") ) y^^^ ( a[") - Ai") ) 



-^(ri)/,(ri)s (n)/,(n) \ (") \ ,,(") ^ \ (") \(")V 



a 7^ 6 = 1,2. 



(n) G (cW,Af)) (6.67) 

for Gn'' and A2^^_]^ due to their continuation for cj"'^'' and A^' respectively. Here we are substituting the 
residue expressions for the eigenvalues given by H6.48|l . 

The nested Bethe ansatz for one and two-particle state presented here contain all information necessary 
to known what happens when a multi-particle state is considered. The sequence of terms in the eigenvalues 
and in the Bethe equations are common for all models and differences will appear only in the last terms. 

7 The multi-particle Bethe state 

Generalization of previous results in order to consider Bethe states with more the two particles follows 
from [2D) where the vector was defined through the following recurrence formula: 



<i>«(Ai, A™) = S(')(Ai) ® <?_i(A2, A„) 

-^-.-(AOE ^ e2(A.) n '4k^ n ^^'h^, - a,)4'^(a.) (t.d 

j=2 VNiNS^'i- ^ '^j) k=2 Xi {Xk - Xj) 



with the initial condition ^q'' = 1, $J'-'(A) — B^''''{X). Here we have used the notation (Xj) to indicate 
the absence of the spectral parameter Xj . 
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In order to proceed with the Bethe ansatz construction we must compute the action of the diagonal 
operators A^'' (A), i = 1, 3 and Tr^ [2?*^') (A)] of the monodromy matrix on the vectors $m . This procedure 
is really very laborious due to the normal-ordered condition but one can do it recursively. Acting with 
A^'^ (A) on (|7.1() we have the following normal-ordered expression 



4')(A)<i>W(Ai,...,A„) = n^''HA,-A)$W(Ai,...,A™)A«(A) 

k=l 

j=i 4 (^j - ^) k=i x\ (Afc - Aj) fc^i 

j=2p=i k=i ^1 '(Afe-Aj) fc^p+i ^(Afe-Aj) 

p-l n{l + i) / \ _ \ A m 

U tt^ , ! n ^('nA.-A,)z(0(A,-A,)A«(A,K^(A,) + ... (7.2) 

= l (Afc Ap) J,^;^ 

where the indices k,k + 1 in the matrix S'|,'^^\ denote the spaces where its action is not trivial. 
Gjp(A, A;, Xj) are matrix valued functions given by 

rWn X X ^ - 4'HA,--A)yff,(A-Ap) g('+i)(Ap-A) ^i^?2(Ap-A,) 
( Aj - A) y]v; ( A - Ap) (Ap - A) y]v, at, ( Ap - A) y]v>, ( Ap - A^ ) 

In (|7.2|) one can easily identify the candidate for the wanted term in the eigenvalue problem and two groups 
of unwanted terms. In each group terms differ by cyclic permutations of rapidities and consequently, by 
the presence of S matrices. 

The action of A'-^ (A) on $m has a similar form 

4')(A)ci>W(Ai,...,A„0 ^ n S ^^U ,W(Ai,...,A„OA<')(A) 

fc=i y]v,iv,(A- Afc) 

j = l VninS^ ^ ^3) k=j + l ^1 l-^i ^ ^feJ fe=l 

™ j-i ^ ^ m S'^'^"'"-' f A ■ — At) 

+i3Ar,_i(A)^5]HS(A,Ap,A,)$|il,(Ap,A,)0pW(Ap)0pW(A,) [] ' 

j=2p=l k=j+ 

m c('+l) \ A i-1 c('+l) ('\ \ \ m 

n'^fc-2,fc-lV'^' ^k) T-r -^fc-Lfcl^p ^fcj -r-r /n , , (n , , , , 

(i+i),. r— 11 ^1+1)77 — 11 ^^'nAp-Afc)z(')(A,-Afc) + ... (7.4) 

fc-=j+i A '(A/ -Afc) fe^p+i '(A, -Afc) j^^i 

fc#pj' 

where we have the following matrix valued functions 

,.(0 M X A vil) l^ \ A AD,. X N f^W 
^ (\ \, \ - - 

3P 



fe=j+i (Aj — Afc) 



. . . y!^!i(A-Ap) yff,(Ap-A,) x^')(A-Ap) yff,(A-A,) 

n,-pl'^:'\^'^i; - -m — - — — ~m — ~ r~ m — ~ — r~~m — 7 — r~ y'-^i 



VnIn, (A - Ap) y'^^] (Ap - Xj) (A - Ap) (A - A^) 

In (|7.4() we have a wanted term and presence of a new group of unwanted terms. 
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Acting with Tr^ [15(A)] on the vector <I>„i the final expression is more cumbersome 



TrjP«(A)]$«(Ai, . . . , A„0 = ci>«(Ai, A„)Tr.[%,i^^ • • ■ %(^^p(0(A)] 

X2 (A — A,„) X2 (A — Ai) 

j=i (A - Ajj k=3+i ^1 (Aj - Afe) 

]=i VninSX - Xj) ^.^-^ x^ (Afe-Aj) j,^;^ 



j=2 p=i fc=i (Afc Apj 



n "^^^ ' n ^^'nA.-A,)z«(A,-A.)4'^(Ap) 



(Aj - Afc) f,=i 



j=2p=i k=i ^1 (Afc — Aj) 
n+i).^ — 11 g^^(Afc - Ap)z^-'(Aj- - Afc) —A\'{\j) + --- (7.6) 



X 



k=J+l ^1 



(Ap - Afc) i^^i Xi (Ap-Aj) 



k^pj 

where we have defined the matrix valued functions 

vWa \ i ^4^^(A-Aj) 4'' (A - Ap) xfjXp - A^) ^ ^ir'^A - Ap) 

Y^p(A,Ap,A,) - [z (A-Ap)^yy- —-^7)- T-^oTT TT^ (/) TT 

4 (A-Aj) 4'(A - Ap) 4'(Ap - Aj) y]v jv,(A-Ap) 



(7.7) 



Here we observe the presence of TZ matrices instead of S matrices as in the expressions for ^i(A) and 
A3 (A). This difference is fundamental for the Bethe ansatz construction. For instance, the trace H7.6|) 
is to be understood as a L + to site inhomogeneous transfer matrix with m inhomogeneous sites coming 
from the Lax operators £^'+^^(A — Afc) = 7^^'+^^ (A — Afc). The solution of this inhomogeneous eigenvalue 
problem is the candidate for the wanted term in the eigenvalue problem of the L site homogeneous transfer 
matrix tl{\). Moreover, the remained terms of ()7.6|l have the exact form to cancel the unwanted terms 
coming from (|7.2(l and H7.4(l . As before, the ellipses used in (|7.2|l . (|7.4|l and (|7.6(l denote normally ordered 
terms containing annihilation operators. 

In our nested Bethe ansatz language the ground {I — 0) for a particular vertex model is prepared by 
a L site homogeneous transfer matrix Ti(A) 



N-2 



tl(A) = Ai(A)+ E ^""(A) + A3(A) (7.8) 

and the multi-particle Bethe state is defined by the linear combination 

4'm(Al, • ■ • , A™) = $™(Ai, \m)Tm |0l) (7.9) 
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where is a vector matrix with (A'' — 2)™ entries 

The corresponding eigenvalue is obtained from the first term on the right hand side of H7.2|l . 17.4(1 and 
(Enjl: 

A,(A|{AJ) = X,iX) n ^(A. - A) + X3(A) [] ' 'f, + X,(A)^^I^^i^^M^ (7.10) 

feii ^j^,yiViv(A-Afc) x«(A) 

Here we have used (|3.9f) and Am^(A, {Ai}|{A|"^'}) is the eigenvalue of the eigenvalue problem for a m-site 
row-to-row inhomogeneous transfer matrix with its Lax operators identified with the matrices TZ^^^ (X—Xt), 
k = 1, ...,TO where Xk are the inhomogenity parameters and {A|^^} are rapidities: 



r^i) {X)T„^ - Tr, (A - A„0 ■ • • ^iV (A - Ai) 

= AW(A,{A.}|{Ar)}).F,„ (7.11) 

where the choice 

.F„ = $a)(AW,...,AW)|oW) (7.12) 

is implicit. 

The remained terms of (17.2(1 . 17.411 and ((7.6(1 multiplied by are known as unwanted terms. 

There are many of these terms but they can be collected in only three different groups. The first group 

A 

contains m terms of the type B{X) (3 $m-i(Aj). To see how is proceeding the cancel in this group, we 
start with j ~ I. The corresponding term has the form 

^^V-^i fe=2 fc=2 a;i^(Ai-Afc) / 

Now we take the limit of t'-^^ (A) at A = Ai in order to identify the product of TZ matrices 
t£) (Ai) = ^hm t£) (A) = ^hm Tr, (A - A„) • • • (A - Ai) 

m— 1 

= n<fe+i(^i-^fc) (7-14) 



fc=i 



Therefore we can use the solution of the second eigenvalue problem ((7.11(1 at A = Ai to see that this term 
will be cancelled provided that 

\0l) = ft '-^^^,Ivi^\Xu{X.}\{X^'}) ft -Jiyj^—-^^ (7.15) 
^2(Ai) ^J-^ z{Xk - Ai) ^J-^ )(Ai _ Xk) 

The remained terms of this group can be written as cyclic permutations of ((7.13(1 . Thus, we can recall 
the relations 1(6. 34(1 - 1(6. 36() to prove that all terms of this group are eliminated provided that 

Xl(A„) ^ (fr z{Xa-xA A^^(Aa,{A,}|{Af^}) 
X2{Xa) yi[z{X,-Xa)j XW{Xa) 

a = l,2,...,m (7.16) 

where 

m 

X^'\Xa) = l[x['\Xa-Xt). (7.17) 
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A 

The second group contains m terms of the type B*{X) (g) <i>m„i(Aj) and they are also canceUed by 
the partial Bethe equations H7.16|l . To accept this statement one can foUow the same steps used in the 

A A 

first group. Finally, the third group contains of the type Bi\i^i{X) ® $m_2(Ap, Xj) are also cancelled by 
1)7.16(1 . Here we would like to stress that the technicalities involving the calculation of the third group of 
unwanted terms are very laborious. 

At this point we have concluded the first step in an algebraic nested Bethe ansatz. The next step 
consists in taking into account the auxiliary eigenvalue problem (|7.11|l . After we have presented the results 
for one and two particle states it is not frivolous to affirm that we will only need to repeat everything 
once more with trivial modifications. It is enough replace L by to and introduce the inhomogeneity 
parameters {A^'"^-*} and the rapidities {A^'-*} for each label I. Indeed this is true until we arrive at the 
last layer I = n — 1, where the models behave differently. 

The last layers for the B^^'' and A'^J^ vertex models are building with the B^^'' Zamolodchivok-Fateev 
model and the A'"^^ Izegin-Korepin model |23[I24| . respectively. For D^^'' models the last layer is the Dj^"* 
vertex model which is mapped in a direct product of two A^^^ six-vertex models. The Cn'^ and A'^J^_^ 

models have their last layer extended to ^ n where we will find the C'i^ and Af ^ six-vertex models, 
respectively. 

From these considerations we can summarize the results in the following way: the eigenvalue of the 
L site homogeneous transfer matrix tl (A) for a m-particle Bethe state is given by 



Al(A|{AJ) = Xi(A)[]z(Afe-A) + X3(A) J] 



X2{X - Xk) 



^J-^ yNNiX - Afe) 
+X2{X) ( ^G(;)(A,{Af-^)}|{Af'}) + r ) (7.18) 



k=l 

/n~2 



\l = l 

where 

G<;>(A, (Af->)iiAf' )) = ^ n ^"'(4" - ^) + ^ n 2 ^'-''il 

^2 l-^J fc=l ^2 'y^) k=iyNiNM~ \ ) 



and we are again using the shorthand notation presented in (|6.31|) and l|t).32|l . i.e., 

m 

X^iX) = [xa{X)]\ Xi'\X}^l[x[!\X-X^t'^) a =1,2. 

m 

X,{X) = [yNN{X)]\ xl\x) = X{y^}}^^^{X-Xt''') (7.20) 



k=l 

The T term in 1)7.18(1 makes the result difference for different models: 

r = G(r^)(A, {A "-)}|{Ar^)}) + n , ^ (^-21) 

for b1^^ and A^^ vertex models, 

T = G(r'^(A,{A("-^^}|{A("-i)}) + G(;')(A,{Af-i)}|{A(")}) (7.22) 



28 



for Cn^ and A'^n_-^ vertex models. For D^"'^'' vertex models the T term is a little bit different due to the 
direct product of two A^^-* six-vertex models: 



T = [X+(A)n^(A^-A)+X+(A)n^(A-A+)] 

m m 

x[Xr(A) n - A) + ^2" (A) n ^(A - Afc )] (7.23) 



k=l k=l 



where 



^(A±) = Xt{X)^Y[xaiX-\r'), (a = 1,2) 

^^^^ > k=i 
,(1) 



{xi,x2} e A\'> (7.24) 
The corresponding sequences of Bethe equations also have a common part 

4'^(A«) " tl^«(Ai'^-Ai'^) ^ ^ ^/i ^ '^ 

/ = l,2,...,n-2. a7^6= l,2,...,m (7.25) 

where we are substituting the residues 

A^^HA«,{Af)}|{Ar)}) = res A^i)(A,{Af)}|{Af+^)}) 

= xr^(o)xr)(A«)nz('+^)(Ar^)-Ai')) 

fe=i 

; = 0,l,...,n-2, a^6= l,2,...,m. (7.26) 

As happened with the eigenvalue sequence (|7.18|) this sequence will end in different ways depending 
on the model: 

^] = TT , n^ (A "''' - Ai"-^^) (7.27) 



for B^^^ and aI?,^ vertex models, 



jT^ilAi^Aj^ ^ fx z{Xt Xf) 
Hx2(A^-A("-)) 



{a;i,X2} e A^'^ (7.28) 

for dI^'' vertex models. For clP and A2^jLimodels we have more two equations: one equation for the 
layer I = n — 1 

4"-'^(Ai"-^)) ,v..("-i)(Ar^)-Ai"-^v H 
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and other equation for the layer I — n 



1 ^ ^ TT ~ ~ f7 29) 

-^3 ) b^a ^1 \\ ~ ^o. ) 2/22 l-^a ~ -^6 ) 

Here we note that the function w(A) used in the B^^^ and A^^n defined in l|().42|l . 

There is a technical point which we already touched in the one-particle case. To solve the eigenvalue 
problem for a L site homogeneous transfer matrix tl{X) with eigenstate ^^(Ai, A2, • • • , A„i) we are left 
with a second eigenvalue problem for a m site inhomogeneous transfer matrix rm'(A, {Xi}) for which the 
vector J-rn must be an eigenstate. However, J>„ defines ^'m as a linear combinations of the components 
of the vector $,„• The dimensions of <i>„i , |0m} and J>„ are suggesting the choice of JF„i as a m-particle 
state in the second eigenvalue problem. The choice of J-m as a r-particle state could give particular 
value for ii r < m and, increasing considerably the number of parameters {a['^} in each layer if 
r > m . However, these different choices are not necessary to fix the rapidities of the states via the Bethe 
equations. 



8 Conclusion 

In this paper the nested Bethe ansatz formulation is used to solve exactly a series of trigonometric vertex 
models based on the non-exceptional Lie algebras. Here a detailed account of this method was described 
in order to complement the results in the literature |17|-|2l)|. 
There are several issues for which this paper could useful: 

a) The off-shell Bethe ansatz - Gaudin theory - Solution of the Kniznik-Zamolodchikov equation. The 
explicit expressions for the eigenvalue problem as presented in this paper define the off-shell Bethe ansatz 
equation. Now, if one can extend the Babujian and Flume formalism 26 for nested Bethe ansatz such 
a issue could be possible. 

b) Graded matrix Bethe ansatz. A recent work jj^l about the vertex models based on superalgebras 
assures the possibility to extend our result for graded vertex models. 

c) A Bethe ansatz with open boundary conditions. The analytical Bethe ansatz for quantum-algebra- 
invariant spin chains |27[I28| gives us the eigenvalues for the models considered in this paper with quantum 
open boundary. The nested Bethe ansatz with diagonal reflection K-matrices was used by de Vega and 
Gonzalez-Ruiz [221 to find eigenvectors and eigenvalues of the An'' vertex models. More recently, the 
nested Bethe ansatz with diagonal K-matrices boundary conditions for the bI^'' vertex model [30] and 

f _____ 

the Ajn vertex models have been considered. Looking at the difficulties of obtaining these results 
we can see how could be useful an unified Bethe ansatz formulation with open boundary conditions. 
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Appendix: ^ vertex models 



For sake of completeness we will describe in this appendix how the aI^"* result is obtained by reduction 
from the our previous results 

The aIP matrix TZ has the form 

Ni Ni Ni 

^^'^ = Xl ^ Eaa «1 Eaa + ^2 E^a «> Epp + X'i^ E^p ® Ep^ 

a^ot' OL^(3 ct<l3 

Ni 



+X4 ^ Eaf3 Ef3a 
a>f3 

where the Boltzmann weights are given by 

xi{X) 
X3(A) = {q'-l) 



a:2(A)=g(e^-l) 



(A.l) 



(A.2) 



Here Ni = n — I + 1. Note that these weights are not labeled by I because they are the same ones for all 
^n~i models. 

After we have identified the Lax operator with this TZ matrix the corresponding monodromy matrix 
can be written as a TV; by Ni matrix 



1 4^ 


i3« 


^2 






^12 




^21 


^22 



The usual reference state in a L site homogeneous lattice is the highest vector of T'-'^ 

CW(A)|Ol)^'^ = 0, i?i'J(A) |0i)(" = 

si"(A)|Oi)('^ ^ {0,|0i)(')}, a^/3 = l,2,...,iV,-l 

where 

xf^(A)-[xi(A)]^ and X^'^(A) = [x2(A)]^ 
Therefore we can write the monodromy matrix as a 2 by 2 matrix 



(0 



D 



B 



(0 

Ni-l 



\ 



^2,Ni-l 



D 



(0 



(A.3) 



a[^\x) i3(')(A) 
C«(A) I?«(A) 



(A.4) 
(A.5) 

(A.6) 



where we identify a scalar yl'^^^(A), two vector ;B(')(A) and C(')(A) with A^; - 1 entries and a A^; - 1 by 
Ni — 1 matrix P'^''(A).The commutation relations among the matrix elements of HA.6(I can be obtained 
using the intertwining relation with 





Xi 











\ 







X4 


X2 










X2 


X3 







V 











SH+i) 


J 



(A.7) 
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where S''-'+^-' is exactly the permutation of (??) with I replaced by Z + 1. In ljA.7|) xi is scalar and 
{x2,X3,X4} are proportional to the identity matrix. 

Now it is easy to see how our previous general results can be reduced in order to obtain the well-known 
results of the An'' vertex models [33]. Removing the third row and the third column of H3.11|l we have 
HA.6p . It means that the entries {Bn,-i, B* , A^^C* ,Cni-i} are vanishing in the cases. Moreover, 
removing from H4.1|) all row and column with entries {Uap} we will get (IA.7|I . Indeed these reductions 
already expected since that the structure of TZ matrices (|2.1|) and (|A.1|I differ (up to normalization) by 
the ya/3(A) terms (|2.3f) . 

Using these reductions in the previous results we are working with the An'' vertex models. For 
instance, one can use the intertwining relation H4.5|l with (|A.fi|l and l|A.7|l in order to derive the following 
commutation relations 



TrjP«(A)]B«(/i) 



X2(fJ- — AJ 



B(')(/i)Tr, 
B'^^\h) ^ B^^\X) 



p(')(A) 



X2{X - IJ) 

^ Xi{X - ^i) 

X2{X - ^) 

50+i)(A-^) 
xi{X - n) 



1(0 

X4:iX - fi) 



B^^^\X)v'^^\p) 



2^2 (A - /i) 

A^^[y,)V^'\X)-A^\x)V^'\n) 



(A.8) 



which are the reductions of (|4.12l) . I5.2|l . (|6.1UI) and (|6.2I) . respectively 

Now we define the normal ordered m-particle vector by the reduction of 17.1|l to 

$«(Ai, A2, . . . , A„0 = 6(')(Ai) ® e(')(A2) B^'HXm) 
Using (|A.8|I one can compute the action of the diagonal elements of (|A.6p on this vector 

4')(A)$W(Ai,A2,...,A™) 

= n ^(^fe - A)$^iHAi, A2, . . . , A™)4"(A) 



(A.9) 



fe=i 



^ ^ ^kl+iiXk - Xj) 



l[z{X,-X,)4\x,) (A. 10) 

fc=i 



and 



Tr,[P«(A)]<i>W(Ai,...,A,„) 

ci>(i)(A„ A^)Tr.[ ^-|'(^-Y^ ■ • ■ ^^Q^^^'^HX)] 

\ \ rii 



^ X2iX- Xj)' 

or we can recall ()7.2|l and (|7.6|l in order to get these expressions. 



Afc) 



(A.ll) 
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The eigenvalue problem on the ground (/ ~ 0) is 

ri(A)*(Ai,...,A,„) = (A(A)+Tr,[P(A)])^'(Ai,...,A„0 

= Al(A|{A,})*(Ai,...,A™) (A.12) 

where the m-particle Bethe state is defined by the linear combination 

*(Ai, . . . , A„) = «>™(Ai, . . . , Xrn)^m \0l) ■ 

Here, is a vector with entries - q-^ — 2, — 1. 

Substituting (|A.10|I and (|A.11|I into l|A.12(l the eigenvalue problem gets the form 



tl(A)^'„(Ai,...,A„) = A:i(A)]^2(Afc-A)*,„(Ai,A2,...,A„) 



fe=i 



where we have made the choice 



+X2(A)^^!^^M^fM^^„(A„ A2, . . . , A™) 

^i(A,) n ^(A'^-Aj) 

-X,(A,) n -(A,-A.)^^5^^^^^^|-^^ (A.13) 

^™ = $W(aW,...,aW)|0™)W (A.14) 

and we are collected the unwanted terms taking into account the cyclic permutation property Ht).36|l . 
Therefore, the eigenvalue is 

A^(A|{AJ) = X,{X) n -(A. - A) + X,(A)^^!^iiMMM^ (A.15) 
k=l ^2 (A) 

provided 

^i(A,-) ^ ( Y\ A^V(A„{AJ|{Af)}) 

where we have substitute Am-* (A^ , {Ai}|{Ap^}) by the residue of Am^(A, {Ai}|{A|^^}) at A = A^ which is 
find by the second eigenvalue problem. The second eigenvalue problem gives us 

A!.'.'(A.(A.}|(Ai"}) _ X!"(A) ;V „ , A!,?(A,(Ai"}|{Af' }) 
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provided that 



= n i"'?v ^(o)n^(Af-Ar) 

k=l,k^j ^V^k > k^j 

where we have substitute A^^'(A^.^\ {A^ ^}|{Af ^}) by the residue of aS)(A, {Af '}|{Af ^}) at A = A^-^^ 
which is given by the third eigenvalue problem, and so on. We follow this procedure till we reach the last 
layer which consists of the six-vertex model whose transfer-matrix diagonalization is well known in the 
literature. Therefore, the eigenvalue of the transfer matrix for the A^,^^ vertex models is given by 



fe=i \;=i k=i 



Ai(A|{AJ) = X,{X) H z{Xk - A) + X,{X) [ ^ G«(A, {Xf-'^}\{X^^}) + l[ziX- A^""^') ] (A.18) 
where 



G«(A, {Af-^)}|{Af)}) = ^ n -(4'^ - ^) (A-19) 

^2 (A) fc=l 



The Bethe equations are 



-^i(Aj) _ _ TT ^(Aj - Afc) (1) 

^ = n |^-(o)n=(Ar'-AS") 

I = l,2,...,n-l (A.20) 

These results complete our study about the nested Bethe ansatz for vertex models based in non-exceptional 
Lie algebras. 
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